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The Collins fragmentation function is extracted from HERMES data on azimutlial single spin 
asymmetries in semi-inclusive deeply inelastic scattering, and BELLE data on azimuthal asymme- 
tries in e^e~-anniliilations. A Gaussian model is assumed for the distribution of transverse parton 
momenta and predictions are used from the chiral quark-soliton model for the transversity distribu- 
tion function. We find that the HERMES and BELLE data yield a consistent picture of the Collins 
fragmentation function which is compatible with COMPASS data and the information previously 
obtained from an analysis of DELPHI data. Estimates for future experiments are made. 

PACS numbers: 13.88.-l-e, 13.85.Ni, 13.60.-r, 13.85.Qk 
O ■ I. INTRODUCTION 

o : 

The chirally odd CoUins fragmentation function describes a possible left-right-asymmetry in the fragmentation 
^ ' of transversely polarized quarks [1-3], the effects of which can be observed in two ways. In e^e~ annihilations it may 
, give rise to specific azimuthal asymmetries [4, 5] where one can access H^'' H^'^ adequately weighted with electric 
' or electro-weak charges. In semi-inclusive deeply inelastic scattering (SIDIS) with a transversely (or longitudinally) 
\ polarized target it may give rise to azimuthal single spin asymmetries (SSA) [6, 7] where it enters in combination with 
' the transversity distribution function h\{x) [8, 9] (or other chirally-odd distribution functions). For these observables 
^ [ transverse parton momenta play a crucial role requiring a careful generalization [10-13] of standard ( "coUinear" ) 
. factorization theorems [14]. 

' The first evidence for the Collins fragmentation function was looked for in a study of SMC data from SIDIS with 
. transversely polarized targets [15] , and in a study of DELPHI data on charged hadron production in e+e" annihilations 
' at the Z°-pole [16]. More recently the HERMES Collaboration has published (and shown further preliminary) data 
\^ ' on the Collins SSA in SIDIS from a proton target giving the first unambiguous evidence that and h'[{x) are 
■ non-zero [17-19], while in the COMPASS experiment the Collins effect from a deuteron target was found compatible 
\ with zero within error bars [20] . Finally, the BELLE collaboration presented data on sizeable azimuthal asymmetries 
in e+e~ annihilations at a center of mass energy of 60MeV below the T-resonancc [21]. 
I ' One question which immediately arises is: Arc all these data from difFcrcnt SIDIS and e+e~ experiments compatible, 

i.e. arc they all indeed due to the same effect, namely the Collins effect? 
^ ■ In order to answer this question one may try to fit the HERMES [18, 19], COMPASS [20] and BELLE [21] data 
'T^ [ simultaneously. While worthwhile trying [22] such a procedure has also disadvantages. When using SIDIS data the 
^ • obtained fit for the Collins fragmentation function is inevitably biased by some model assumption on the transversity 
distribution, and the — presently not fully understood — scale dependence of is neglected. 

For these reasons we shall follow a different strategy and extract separately from HERMES [19] and BELLE 
' [21] data, and compare then the results to each other and to other experiments. 

By focussing on adequately defined ratios of to Di ( "analyzing powers" ) we observe that the indications from 
DELPHI [16] and the data from HERMES [17-19], COMPASS [20] and BELLE [21] are compatible with each other. 
Experience with spin observables where radiative and other corrections are known, see e.g. [23, 24], suggests that such 
analyzing powers could be less scale-dependent than the respective absolute cross sections. 

For the analysis of the HERMES data we use [25] predictions for h1{x) from the chiral quark-soliton model [26] and 
assume a Gaussian distribution of transverse parton momenta. Our results are in agreement with the first extraction 
of from the HERMES data [27], where different models were assumed.^ On the basis of a comparison to Ref. [27] 
we discuss which of the observations made here and in [27] are robust, and which are likely to be model-dependent. 

We shall furthermore present some estimates for future experiments in SIDIS and at e~^e~ colliders. The comparison 
of these estimates to the outcome of the experiments will help to solidify our understanding of the Collins fragmentation 
function and the transversity distribution. 

The SIDIS data from HERMES and COMPASS [17-20] also provide the first unambiguous evidence for the equally 
interesting Sivers effect [38-41]. We refer the reader to Refs. [27, 42-47] for studies of this effect. 



X 



^ Previously, the preliminary result [15] and longitudinal target SSA [28-31] data were used to extract from SIDIS [32, 33]. However, 
in these studies twist-3 effects [34] were not or only partially considered. Still earlier attempts to learn about the Collins effect from 
SSA in the hadronic processes pp^ or ppl — > nX [35] were made in [36] , but recent studios indicate that the Collins effect is not likely 
to be the dominant source of SSA in these processes [37] as it was assumed in [36]. 
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II. COLLINS EFFECT IN SIDIS 



For reasons which will become clear in the sequel we start our study with the Collins effect in SIDIS — although 
this case is more involved because here in addition to also the unknown transversity distribution h'^{x) enters. For 
the latter we use here predictions from the chiral quark-soliton model [26] which provides a consistent and successful 
[48; 49] field theoretic description of the nucleon. The distribution functions /f (x) and 5° (x) computed in this model 
at a low normalization point [50] agree to within (10-30)% with parameterizations [51]. One may assume the same 
accuracy for the model predictions for transversity. which arc inbetwcen the popular working assumptions, namely 
that at a low scale h1{x) saturates the SofFcr bound |ft.i(a;)| < ^{ff + gi){x) [52] or is equal to gf{x). 

The way the HERMES [18, 19] (and BELLE [21]) data were analyzed requires to assume also a model for transverse 
parton momenta. We assume the distributions of transverse parton momenta in the transversity distribution and 
transverse hadron momenta in the Collins fragmentation function to be Gaussian, and take the respective widths 
(p^J and (K|fJ to be flavour and x- or z-independent 



h1{x,p^T) 

ff^(^,K|) 



hUx) 



exp(-p2,/(p2j) 



exp(-K2,/(K2,J) 



(1) 



The Gauss Ansatz can be considered at best as a crude approximation to the true distribution of transverse momenta 
[53]. For example, it docs not yield a 1/Phi_ power-like suppression of the Collins SSA expected at large transverse 
hadron momenta [2]. However, this Ansatz gives a satisfactory effective description in several hard reactions, provided 
the transverse momenta are much smaller than the hard scale of the process [54] . This is the situation in the HERMES 
experiment where the mean transverse hadron momenta {Ph±) ^ 0.4 GeV ^ yiQ^ ~' 1.5 GcV. 

In fact, the HERMES data on {Ph±) can be described satisfactorily 
assuming for fi{x,p'^) and Di{z,K.'^) Ansatze analog to (1) with [45] 



(p^J = 0.33GeV^ and (K^J = 0.16 GeV^ 



HADRON PRODUCTION PLANE 



(2) 



Noteworthy these numbers are in good agreement with results [43] inferred 
from EMC data [55] on the so-called Cahn-effcct [56] . All that can be said 
about the Gaussian widths of the transversity distribution and Collins 
bmction is that positivity conditions require [57] 
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(3) 



which are necessary conditions (but not sufficient for H^, see [45] for the 
discussion of the analog case of the Sivcrs function) . 

In the Gaussian model (1) the Collins SSA as measured by HERMES 
[18, 19] is given by (we neglect throughout the soft factors [11, 12]) 



FIG. 1: Kinematics of the SIDIS process 
Ip I'hX and the definitions of azimuthal 
angles in the lab frame. The target polar- 
ization vector is transverse with respect to 
the beam. 



_L(l/2)o 



(4) 



where x = QV(2P • g), z = P • Ph/P ■ q, 



-q^ and P is the target momentum. Other momenta and azimuthal 



angles are defined in Fig. 1. The (l/2)-moment of the Collins function is defined as and satisfies the inequality 



H 



±{l/2)a 



^Ht^{z,KT)<^D1{z) 



2zm^ 2 

which follows from (3), see [57]. In Eq. (4) the dependence on the Gaussian model is contained in the factor 

1 



Gauss 



l + z2 {pl)/{Kl^) 



(5) 



(6) 



Within the Gaussian model Eq. (4) is just one convenient way of writing the result. Equivalently one could write 



jj „±(l/2)a. X 

iiGauss-ni [Z) 



^Gauss-^1 



.(l)a 



(z) , aca 



{^'hi)/z' 



(7) 
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Here the (l)-transverse moment of the Colhns function is defined as in Eq. (5) but with an additional power of 
|Kt|/(2to^) in the weight, and acauss is defined in complete analogy to the notation of Ref. [45]. 

Had the Collins SSA been analyzed as {s'm{(l) + (j)s)Ph±/ {zniTr)) , i.e. with an additional power of transverse hadron 
momentum Ph± in the weight, the result would be given by an expression analog to Eq. (4) but with i?Gauss 1 
-^(^/^^ replaced by H^^^^ independently of any model of transverse parton momenta [7]. It was argued that 



and H 

adequately weighted SSA might be less sensitive to Sudakov suppression effects [58] 
analyzed in this way were presented in [17]. 



Preliminary HERMES data 



By defining A^^ = A 



sin{4>+4>s) • 



Y,a el xfl Dl'"" and using HERMES data [18, 19] on 7r± we can rewrite (4) as 



C/T,7r 



hxhf - 



/25, 



Gauss 



H 



xhf J VS-Bcauss-ff: 



±(l/2)fav 
i(l/2)unf / ' 



which — with our chosen model for h1{x) — can be inverted to give unambiguous results for Bq 
favoured and unfavoured Collins fragmentation functions are defined as, schematically: 



auss-^l 



±(l/2)a 



u/tt 



Here — in the study of SIDIS data — we neglect the effects of strange 
and heavier flavours, which is justified at the present stage because the 
corresponding distribution functions are rather small. 

As a first step we extract in this way from the HERMES prelim- 
inary data [19] on the ^-dependence of the Collins SSA the quan- 
tity (2_BGauss-f^i'''^''^''") which is averaged over z within the HERMES 
cuts 0.2 < z < 0.7. Using for the chiral quark soliton model 

[26] and for /"(x) and D1{z) the LO-parameterizations [51, 59] at 
(g2) = 2.5 GeV^ which is the average scale in the HERMES experi- 
ment, we obtain the results shown in Fig. 2. As the SSA for 7r+ and 
7r~ are independent observables the statistical errors shown in Fig. 2 
are not correlated. 

The (2i3Gauss-ffi"'^^^'''^) are expected to be a;-independent — pro- 
vided the used models, the Gaussian Ansatz (1) and the model [26] 

for h1{x), are appropriate. In fact, the extracted (2i3Gauss-ffi''^^^'°) 
are compatible with this expectation and can be fitted respectively to 
constants with reasonable per degree of freedom (x^of) 



TT± (1/2) favv 



(8) 



The 
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(a) 
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(b) 



(2i?Gauss^^^l 
(2i?Gauss-ffi 



_L(l/2)fav 



_L(l/2)unf 



)= (3.5 ±0.8)%, xLf = 0.6, (10) 
) =-(3.8 ±0.7)%, Xdof = 1.3. (11) 
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FIG. 2: The quantity 2Bgs.ussH^^^^'^'''" aver- 
aged over z, i.e. practically the weight of hi{x) 



in the Collins SSA .4" 



{x) in Eq. (4), vs. 
X as extracted from the preliminary HERMES 
data [19]. This quantity does not show any 
significant x-dependence — as expected, see 
text. 



These constants are indicated in Fig. 2. The published HERMES data 
which have a lower statistics [18] yield similar values but with larger 
uncertainties and worse xlof ■ 

Let us remark that at small x a deviation from the expected be- 
haviour (2i?Gauss^^i"^^^^''"(a;)) = constant would not be surprizing for 
two reasons. First, the description of distribution functions in the 
small- a; region x < 0.05 is beyond the range of applicability of the 
chiral quark-soliton model [50]. Second, in the lowest x-h'm of the 
HERMES experiment (Q^) = 1.3 GeV [31] could be at the edge of the 
applicability of the factorization approach. 

However, the lack of a noticeable a;-dependence of (2i3Gauss^fi^'^^^''") 
indicated in Fig. 2 gives certain confidence that the model [26] for hi{x) 

used here is — considering the accuracy of the data [18, 19] — reasonable. The results in Eqs. (10, 11) are, of course, 
model-dependent but Fig. 2 indicates that the systematic uncertainty due to model-dependence is less dominant than 
the statistical error in (10, 11). 

We also remark that the results in Fig. 2 and Eqs. (10, 11) are not specific to the Gaussian Ansatz: Any analysis 
of the data [19] in which the factorized Ansatz hi{x, pf.) = /i°(x) G{p'^) is assumed with G(p^) some function of 
(and h^ix) is taken from [26]) would yield the same result (10, 11). (Then, however, this z-averaged quantity would 
be given by some different theoretical expressions.) 

The results (10, 11) are unexpected from the point of view of studies [33, 60-63] of longitudinal SSA [28, 29]. In 
order to understand these data it was sufficient to assume favoured flavour fragmentation only and to neglect ifj^""^ 
completely, see [63] for a review. Later it became clear that these SSA are dominated by subleading twist effects 
[31], and are theoretically more difficult to describe [34]. Moreover, we now learn that the unfavoured fragmentation 
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^sin(Ms)(^j for proton (a) Agil-'^+^sjix) for proton (b) Aj^*t:(x) for deuteron (c) Ag!f*<:(x) fordeuteron (d) 
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FIG. 3; The Collins SSA A''^"^'^^''"'\x) as function of x. The preliminary HERMES data are from [19], the COMPASS data 
are from [20]. The theoretical curves are based on the fit in Eqs. (10, 11) and predictions for the transversity distribution from 
the chiral quark-soliton model [26]. Notice that different sign conventions are used in [19, 20]: A!'^^'^ {x) = —A^^'^^'^^\x). 



in the Collins function cannot be neglected. Instead, in order to explain the HERMES data [18, 19] the favoured 
and unfavoured Collins fragmentation functions must be of similar magnitude and opposite sign. The string frag- 
mentation picture provides a qualitative understanding of this behaviour [64] as do quark-hadron-duality motivated 
phenomenological considerations [27]. 

The different behaviour of the Collins function compared to the unfavoured fragmentation function becomes more 
evident by considering the analyzing powers defined as 



auss 



_L(l/2)fav, 



fav\ 



(7.2 ± 1.7)% 



(2-BGaussffi 



±(l/2)unf. 



HERMES 



-(14.2 ±2.7)% 



(12) 



HERMES 



The positivity bound (5) and the fact that iJcauss < 1, see Eq. (6), require the absolute values of the above numbers 
to be smaller than unity which is the case. Thus, the extracted results satisfy positivity. 

Figs. 3a and b shows how the fit (10, 11) describes the preliminary HERMES data [19] on the Collins SSA from 
the proton target. As can be seen from Figs. 3c and d this fit is also in agreement with the COMPASS data which 
show a compatible with zero Collins effect from a deuteron target [20] . 

In principle, one could extract in an analog way also the z-dependence of the product Bcauss-ffi^^^^^^'^- However, we 
refrain from doing this here for the following reasons. First, due to the appearance of the unknown Gaussian widths 
in i?Gauss, see Eq. (6), we could not conclude anything on the z-dependence of H^^^^'^'^°' . Second, the study of the 
z-dependence requires to average over x within the HERMES cuts 0.023 < x < 0.4 which includes small x where the 
model [26] is not applicable. We postpone the study of the z-dependence of the Collins SSA for Sec. IV. 



III. COLLINS EFFECT IN e+e ANNIHILATION AT BELLE 



In the process e+e^ /11/12X where the two observed hadrons belong 
to opposite jets the Collins effect gives rise to a specific azimuthal dis- 
tribution of hadron hi along the axis defined by hadron /12 in the e^e" 
center of mass frame [4, 5], sec Fig. 4. Assuming the Gaussian model this 
azimuthal distribution is given by 



d^g(e+e-^ hih2X) _ 3a^ 
dcos02dzidz2d0i 

+ COs(20i)sin2 02CGaus.E^'^l^*'^''"(^l)^l^^'^'^'(^2) 



(13) 



where is the center of mass energy square of the Icpton pair and 
Zi = 2Efi./Q while Ccauss is defined as 



Cg auss 



16 Z1Z2 



(14) 




FIG. 4: Kinematics of the process e"'"e~ 
hih2X and the definitions of azimuthal an- 
gles in the e"'"e~ rest frame. 



In principle, one can rewrite the result in (13) within the Gauss model also in terms of or H^^^'' with some 
different functions instead of Ccauss- 

Experimentally it is convenient to normalize the expression (13) with respect to its (/)i-average and to define the 
following observable [21] 



Ai = l+ cos(2(/)i) 



(sin 



rUl/2)a 



(^2) 



(1 + C0S2 6*2) 



T.a^lDl{zi)Dl{z2) 



(15) 
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Here the average over 02 is understood in the range of acceptance of the BELLE detector [21]. 

Notice that hard gluon radiation also gives rise to a cos(2(^i)-depcndcncc as do detector dependent effects. These 
effects, however, are flavour independent and — as long as the coefficient of the cos(2(^i)-niodulation is not large — 
one can get rid of them in the following way [21]. One considers the double ratio of , where both hadrons are 
pions of unlike sign, to A^, where /11/12 are pions of like sign, i.e. 



A^ 

= l + cos(20i)Pi(zi,Z2) . (16) 



The observable Pi{zi, 22) is given by 



(sin^ 02) 



(1 + COS^ 02) 



r Tj-L(l/2)fav, X jj_L(l/2)fav. x , jj-±(l/2)unf . x jj_L(l/2)unf / x 

Si?! (zi)H-^ {z2) + 7H^ {zi)H^ (Z2) 

bDf^{zi)Df^{z2) + 7L»5'"f(zi)i:'™f(z2) 



3^^^(l/2)fav(^^^^±(l/2)unf (^^^ + bUt^' f'^''''\z,)Ht^"^^"^^ {Z2) + 2g,^<^/^^""^(zi)g,^<^/^'""^(z2) 

bDf^{zi)DY'^{z2) + bDY'^{zi)Df^{z2) + 2Df'^{zi)DY'^{z2) 



(17) 



up to higher order terms in the Collins function. The systematic error of the double ratio method was estimated to 
be smah [21]. 

In Eq. (17) it is assumed that the Collins fragmentation of s- and s-quarks into pions is equal to the unfavoured 
fragmentation function defined in Eq. (9). At the high energies of the BELLE experiment just below the threshold 
of 5-quark production this is a reasonable and commonly used assumption for [65] . We assume it here to be valid 
also for the Collins function. [D^^^ and Z?""^ in (17) are defined analogously to Eq. (9).) Charm contribution does 
not need to be considered in Eq. (17) since the BELLE data are corrected for it [21]. 

In order to obtain a fit to the BELLE data we adopt the LO-parametcrization [65] for -D" (z) at ~ (10.52 GcV)^ 
and choose the following simple Ansatz 

Ht'^^'^^\z)^CazDl{z) . (18) 

The two free parameters Ctav and Cunf introduced in (18) can be well fitted to the BELLE data [21]. We explored 
also other Ansatze proportional, for example, to z^D°(z), D'l{z){l — z) or zD\[z)(\ — z) but none of them gave 
satisfactory solutions. The best fit has a Xdof ~ 0-6 and is demonstrated in Fig. 5a in the Ctav-Cunf-plane. Two 
different, equivalent, best fit solutions exist. The reason for this is that the expression for Pi{zi,Z2) in Eq. (17) is 
symmetric with respect to the exchange Cfav Cunf in our Ansatz, and manifests itself in Fig. 5a where the two 
solutions are mirror images of each other with respect to the axis defined by Cfav = Cunf- 

What can unambiguously be concluded from Fig. 5a is that the BELLE data require the Collins favoured and 
unfavoured fragmentation functions to have opposite sign — as in the HERMES experiment. On the basis our study 
of the Collins effect in SIDIS we are tempted to select the solution with Cfav > and Cunf < in Fig. 5a as the 
appropriate one. Our result is thus 

Cfav = 0.15, Cunf- -0.45 (19) 

and the resulting best fits and their 1-a regions are shown in Fig. 5b for z > 0.2 which is the low-z cut in the BELLE 
experiment. The results satisfy the positivity condition (5). Notice that the errors of the favoured and unfavoured 
Collins functions are correlated. 

In Fig. 6 the BELLE data [21] are compared to the theoretical result for Pi(zi,Z2) obtained on the basis of the best 
fit shown in Fig. 5b. Here the error-correlation of the fit results for the favoured and unfavoured Collins function is 
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Hj-*'"'(Z) at BELLE (b) 
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FIG. 5: a. The two best fit solutions for the parameters d in the Ansatz (18) (indicated as discrete points) and their respective 
l-CT regions as obtained from a fit to the BELLE data [21]. The solutions are symmetric with respect to the fine Cfav = Cunf 
indicated by a dashed fine. b. The best fit for H^'^^^'^^"' {£) resulting from Fig. 5 with the choice H^^^'' > and ffj^""* < as 
suggested by the analysis of the HERMES experiment, see Sec. II. 
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P(Zj, Zj) tor 0.2 < Zj < 0.3 (a) P(Zi, * ^1 (b) P(Zl! "-^ ^1 ^ (C) P(2l> f"' 0.7 < Zj < 1 (d) 




0.2 0.4 0.6 0.8 Zj 0.2 0.4 0.6 0.8 Zj 0.2 0.4 0.6 0.8 Z2 0.2 0.4 0.6 0.8 Zj 



FIG. 6: The observable Pi{zi, Z2) as defined in Eqs. (16, 17) for fixed zi-bins as function of Z2. The data are from the BELLE 
experiment [21]. The theoretical curves are obtained on the basis of the fit result shown in Fig. 5. 



taken into account and the resulting l-cr error band is more narrow than in Fig. 5b. The description of the BELLE 
data is satisfactory — as can be seen from Fig. 6. (Notice that Pi(zi, Z2) is symmetric with respect to the exchange 
Zi Z2 such that the bins with zi ^ zi can be combined — as was done in the BELLE analysis [21] and in our 
fit procedure. Here, for a better overview, these bins are presented separately — whereby we disregard that strictly 
speaking the statistical error bars for bins with z\ ^ Z2 should then be multiplied by ^/2.) 



IV. ARE BELLE AND HERMES DATA COMPATIBLE? 



In order to compare the Collins effect in SIDIS at HERMES [18, 19] and in e+e--annihilation at BELLE[21] it is, 
strictly speaking, necessary to take into account the evolution properties oi . However, in order to get a first rough 
idea — which is sufficient at the present stage — we instead consider ratios of to _DJ, which may be expected to 
be less scale dependent. For example, integrating the BELLE fit result in Fig. 5b over the range of HERMES z-cuts 
0.2 < z < 0.7, we obtain the following analyzing powers 



_(l/2)favv 



fav\ 



(5.3 . . . 20.4)% 



(2H- 



_L(l/2)unfv 



BELLE 



(3.7 ... 41.4)% 



(20) 



BELLE 



Comparing the above numbers (again the errors are correlated) to the result in Eq. (12) we see that the effects at 
HERMES and at BELLE — as quantified in Eqs. (12. 20) — are of comparable magnitude. The central values of the 
BELLE analyzing powers seem to be systematically larger than the HERMES ones. This could partly be attributed 
to evolution effects. However, notice that in the HERMES result (12) in addition the factor iJcauss < 1 enters, which 
tends to decrease the result. Thus, the HERMES [18, 19] and BELLE [21] data seem in good agreement. 

Encouraged by this observation let us see whether we can describe the HERMES data on A^!^^'^+'^''^(z) on the basis 
of the z-dependence of concluded from the BELLE data [21]. For that let us assume a weak scale-dependence 
not only for z-averaged ratios — as we did above — but also 



.(l/2)a 



(^) 



D-[z) 



±(l/2)a 



(^) 



BELLE scale 



D-{z) 



(21) 



HERMES scale 



Nothing is known about the Gaussian widths of the transversity distribution and the Collins function which enter the 
factor -Bcauss in Eq. (6). Let us therefore assume their ratio to be similar — let us say to within a factor of two — 
to the corresponding ratio of the Gaussian widths of the unpolarized /{'(x,p|n) and D°(2,K|n) in Eq. (2), i.e. 



1 < 



<4 



(22) 



It is gratifying to observe that within the range (22) the factor Bcauss varies moderately between 10% (at small 
z ^ 0.3) and 25% (at large z ~ 0.6). Taking into account the l-cr uncertainty of the BELLE fit shown in Fig. 5 and 
the variation in (22) we obtain the result in Fig. 7. The description of the preliminary HERMES data [19] on the 
z-dependence of the Collins SSA obtained in this way is satisfactory. 

The good agreement observed in Fig. 7 gives further support to our observation that the HERMES [18, 19] and 
BELLE [21] data are in good agreement. Furthermore, we are lead to the conclusion that the assumption of weak 
scale-dependence (21) is reasonable — given the accuracy of the data [18, 19, 21]. 

Finally, we recall that in the expression for J^^'^'^'^^\z) certain integrals over x enter which extend down to 
low X = 0.023 at HERMES [18, 19] where the used predictions for h1{x) from the model [26] are at the edge of 
(if not beyond) their range of applicability. This prevented us from making a direct quantitative extraction of the 
z-dependence of the analyzing power from SIDIS data in Sect. II. However, here we used the model results [26] 
merely for sake of a qualitative comparison. The good agreement in Fig. 7 indicates that the uncertainty due to 
model-dependence is within the statistical error bars of data and uncertainties of the Gaussian Ansatz as estimated 
in (22). 
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V. ARE BELLE DATA AND DELPHI RESULT COMPATIBLE? 



Another interesting question is whether the BELLE and HERMES data are compatible with the rcsuhs obtained 
from an analysis of DELPHI data [16] where azimuthal asymmetries in e^e~ — *■ hh'X at the Zp-peak were studied. 
At DELPHI h, h' were leading charged hadrons in opposite jets. Also in this process the Collins effect gives rise to 
a cos(2(/)i)-asymmetry analog to (13) but with electric charges replaced by electro-weak ones and with modifications 
due to the finite mass and width of Zq [4, 5]. The cross-section differential in (j>i can be written as [4, 5] 

daie+e" ^ hh'X) ^( , ^ „ \ 

^ — = Po 1 + cos(20i) P2 (23) 



with 



' (l+COS^e,) /dzi/dz2 E.,vEa,.C?<''(zO<'''(z2) 

The c° are defined as (the "+" sign refers to i = 1, the "— " sign to i = 2) 



c« 2 = (r| - 2e« sin^ Ow)' ± [T^f (25) 

with the wcak-isospin ^T§^^, = r| = r| = -i. 

The result obtained from the analysis of the DELPHI data reads [16] 

P2,DELPHI - -(0.26 ± 0.18)% . (26) 

In order to have an idea whether the result (26) is compatible with BELLE let us make the following rough estimate. 
We assume a weak-scale dependence similarly to (21) and take the charged hadrons observed at DELPHI to be pions 
and kaons (which were, in fact, the most prolific particles). Then we obtain from the best fit to BELLE data with 
the unpolarized fragmentation functions at = M"^ from [65] 

-P2. from BELLE «-( 0^ + OWZ + QMA + 0T03 + 0^ 0i)24 )%:=- (0.06 ... 0.29)% . (27) 

u,d.s — ^TT.TT u.d,s — >TT,K u.d,s — >K,K c^b — *7r, -tt c,b — >7Z,K c^b — ^K,K 

The different contributions are explained as follows. The contribution "u, d, s ^ tt, tt" is due to Collins fragmentation 
of light quarks into pions. This is the only contribution which is determined by the best fit to the BELLE data without 
additional assumptions. Since pions and kaons are considered as Goldstone bosons of spontaneous chiral symmetry 
breaking one is lead to the conclusion that in the chiral limit [61] 

hm Ht'-"^^'''^{z)^ hm i/i^(^/^)"/^(z) . (28) 

Notice that the analog relations with Hi or H^^^"^ would contain explicit factors of hadron masses. In the real world 
the unpolarized fragmentation functions for pions and kaons are found to be far off the analog chiral limit relation. 
Thus, one may expect the real- world Collins functions to be similarly far off the chiral limit relation (28). However, if 
the way off the chiral limit to the real world proceeded in a spin-independent way, one still could expect the relations 
(18, 19) between the Collins and unpolarized fragmentation functions for pions to be valid approximately also for 
kaons. This is what we did for the sake of our rough estimate. Future SIDIS and e'^e~ data on kaon asymmetries 
will reveal to which extent these assumptions arc justified.^ Finally, to estimate the contribution due to the Collins 



^sin(i»+?s)^^l for proton (a) Aj5.'***s'(z) tor proton (b) 




0.2 0.3 0.4 0.5 0.6 z 0.2 0.3 0.4 0.5 0.6 i 



FIG. 7: The Collins SSA Ap^^'''~^'''^\z) as function of z. The preliminary HERMES data are from [19]. The theoretical curves 
are based on the fit of Hi to the BELLE data shown in Fig. 5b under the assumption (21). For hi the model prediction [26] 
is used. The dashed lines indicate the sensitivity of the SSA to the unknown ratio of the Gaussian widths of hi and in the 
range (22). The description of the preliminary HERMES data [19] is satisfactory. 



^ The relations (18, 19) cannot hold for all hadrons. Otherwise it would be impossible to satisfy the Schafer-Teryaev sum rule [73]. 
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fragmentation of heavier quarks into pions (and kaons) we assume the relations (18) (and their kaon-analoga) to hold 
also for charm and bottom with Cc = Cb = Cunf from (19). These latter assumptions seem natural — at the Z'^-peak, 
where differences in the fragmentation of non-valence quark-flavours q into a specific hadron might be expected to 
be of 0{mq/Mz)- In fact, one finds for pions Z?""^ = Df « D'^^'^ > D\^^ at = M| and qualitatively similar for 
kaons [65]. But the above assumptions need not to be correct and, for example, the Collins fragmentation of heavier 
flavours could be simply zero. This gives rise to the uncertainty in the final result in Eq. (27). 

The estimate (27) — even if we consider only the Collins effect of light quarks — is in reasonable agreement with 
the result from DELPHI in Eq. (26). We refrain from quoting in Eq. (27) the error to due to the l-cr uncertainty of 
the BELLE fit as our result has presumably even larger theoretical uncertainties. Also it is important to keep in mind 
that the DELPHI result (26) has unestimated (and presumably sizeable) systematic uncertainties. What is important 
at this point, however, is that both numbers (26, 27) are of comparable order of magnitude — indicating that the 
DELPHI result [16] and the BELLE data [21] could be well due to the same effect. 



VI. REMARK ON THE PRELIMINARY SMC RESULT 



A preliminary resuh^ for the CoUins SSA An = was ob- 

tained for charged hadrons from an analysis of SMC data on SIDIS from 
transversely polarized proton and deuteron targets at (Q^) ^ 5 GeV^ and 
(x) - 0.08 with (z) - 0.45 and {Ph±) ~ (0.5 - 0.8) GeV [15]. A non-zero 
Collins SSA for positive hadrons from a proton target was reported, the 
other SSA were found compatible with zero within error bars, see Fig. 8. 

The charged hadrons at SMC were mainly tt^ [15]. Neglecting the con- 
tributions from other hadrons we obtain the results shown as open circles 
in Fig. 8. Our results are not in conflict with the preliminary SMC results 
except for the only case where indications for a non-zero SSA were seen, 
namely for tt"'" from a proton target. Is this a discrepancy and if so, how 
could one explain it? 

Could, for example, the Collins effect of hadrons other than pions be 
of importance in the SMC kinematics? On the basis of Schafer-Teryaev 
sum rule [73] one may, in fact, expect cancellations between contributions 
of different hadrons in SIDIS where the Collins effect is linear in H^. 
(In e"'"e~-annihilations the effect is proportional to H^'^H^^ and different 
signs of the Collins functions for different hadrons cancel.) 

It is difficult to address these questions — in particular because the results 
and their systematic uncertainties are not known. Future data on the Collins 
situation. 



20 



Ajy in % 



10 



-10 



-20 



'- c 


SMC preliminary ■ 
our estimates G 


: o 1 


h h 



proton deutron 



FIG. 8: The Collins SSA Ajv = 
— '^"^'^^'^ in the production of charged 
hadrons 



in SIDIS from different targets. 
The squares denote the preliminary SMC 
data [15] . The open circles show our rough 
estimates. 

of Ref. [15] retain a preliminary status 
effect of other hadrons may clarify the 



VII. COMPARISON TO THE EXTRACTION OF REF. [27] 

In Ref. [27] Ht''^'^\z) was extracted from the HERMES data [19] in a similar way as in Sec. II but with some 
different model assumptions. It is worth to inspect these differences in some more detail. This may shed some light 
on the model dependence of the two approaches. (Notice the different notation 5(p-^'^\z) = —2H^^^^^\z) in [27].) 

There are two main differences between the approach of [27] and ours. First, in [27] it was assumed that the 
transverse momentum of the hadrons produced in the SIDIS process is solely due to the fragmentation process, i.e. 
the possibility was disregarded that partons are emitted from the target with non-zero intrinsic transverse momenta. 
This corresponds to setting h1{x,p^) = h1{x) i5*^^-'(pt) and is equivalent to our Gaussian Ansatz (1) in the limiting 
case p^^ ^ in which Bcauss ^ 1 in Eq. (4). Second, for the transversity distribution the saturation of the Soffcr 
bound [52] at a low input scale was assumed as discussed in [72]. 

The extracted analysing power in SIDIS is valid in any model approach where the factorized Ansatz /i°(x,p|n) = 
hi{x) G{pt) is assumed (see Sec. II). This was done in [27] and therefore, it is possible to compare our results (12) 
directly to the fits of Ref. [27] which yield 

r (6.1 ± 0.8)% for set I 
\ (6.1 ± 0.4)% for set H, 

J -(11.0 ±1.3)% for set I 
\ -(11.5 ±1.4)% for set H. 



fav\ 



{2H, 



_L(l/2)unf, 



Ref. [27] 
Ref. [27] 



^ Notice the different convention used in Ref. [15] where the azimuthal distribution of hadrons was analyzed as function of the "Collins 
angle" ipc = fph + <j>S ~ f - Hence, the opposite sign. The same sign convention is used at COMPASS, see caption of Fig. 3. 
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Wc stress that (12) and (29) show the same quantity but in different models of transverse hadron momenta. Sets I, 
II refer to different Ansatze in [27] for the CoUins fragmentation function. We observe a good agreement with (12). 

We remark that the Ansatze of [27] are of the type H^^^^^\z) cx z{l — z)Di{z). The factor (1 — z) in this Ansatz 
is based on general considerations in QCD [2] and is "needed" in the approach of [27] to reproduce the tendency of 
the Collins SSA to decrease (or, at least, not to grow) with increasing z at HERMES [18, 19]. Interestingly we find 
that such an Ansatz is not suited to describe the BELLE [21] data which extend to higher z than the HERMES 
data [18, 19]. In our Gaussian model for transverse parton momenta, where the mean transverse momentum of the 
produced hadrons is taken ^-independent, the "needed" decrease of the Collins SSA at larger z is provided by the 
factor i?Gauss (6). Thus, although both approaches describe the z-dependence of the Collins SSA comparably well, the 
Collins fragmentation functions extracted here and in [27] differ in particular at large-z. This point nicely illustrates 
the model dependence both of our results and those of Ref. [27]. Both approaches are not meant to give a fully realistic 
account of transverse parton momentum dependent process — but are suggested as effective descriptions which are 
useful at the present stage and will have to be improved upon the impact of future, more precise data. 

Notice that the good agreement of the analyzing powers (12) and (29) indicates that there is little sensitivity to 
the choice of the model for the transversity distribution function. We shall come back to this point below in Sec. IX. 

VIII. COMPARISON TO MODEL CALCULATIONS OF THE COLLINS FUNCTION 

The appearance of any SSA requires the presence of interfering amplitudes with different relative phases [66]. The 
Collins function is losely speaking the imaginary part of such interfering amplitudes and this was explored in model 
calculations in Ref. [2] and further in [67-71]. The common spirit of these calculations is to generate the required 
imaginary part by a one- loop diagram — using effective quark degrees of freedom, and pions or (abelian) gluons. 
Such perturbative calculations do not aim at providing a realistic model for the dynamics of the non-perturbative 
fragmentation process. Rather their importance is to indicate the existence of the mechanism [74]. It is therefore 
not surprizing to observe that the models [67-71] have difficulties to reproduce the quantitative properties and in 
particular the flavour dependence of H^, though this could possibly be achieved by exploring the freedom in the 
choice of model parameters. Further studies in this direction are required. 

Presumably better suited to catch features of the non-perturbative fragmentation process is the approach of [64] 
based on the string fragmentation picture. Though quantitatively plagued by uncertainties in the choice and tuning 
of parameters this approach a priori predicted the favoured and unfavoured Collins fragmentation functions to be of 
similar magnitude and opposite sign. Explanations of this observation have been given in [27] on the basis quark 
hadron duality considerations. 

IX. HOW MUCH DO WE KNOW ABOUT THE TRANSVERSITY DISTRIBUTION? 

After the discussion of the picture of the Collins function emerg- 
ing on the basis of experimental information [16-21], it is worthwhile 
asking how much do we really know about the transversity distribu- 
tion. In other words, to what extent do HERMES and COMPASS 
data [17-20] on the Collins SSA in SIDIS constrain h1{x)7 

In order to gain some insight into this question let us make the 
following exercise. We repeat the extraction of the analyzing power 
from the preliminary HERMES data [19] and use for /i" (x) the pre- 
diction from the chiral quark soliton model [26] but for hf{x) we use 
instead the model 

/^^(x)=r;^^(x)|p^f [26] • (30) 

The extracted result for (2_BGauss-ffi'") is then a function of the pa- 
rameter r and its sensitivity to r is shown in Fig. 9. Our actual 
result based on the model [26] is given for r = 1 and its 1-a regions 
are indicated in Fig. 9 by dashed lines. We make the remarkable ob- 
servation that the extracted results for the favoured and unfavoured 
{2BqiiussH^) are nearly insensitive to r. Varying |ft,f(.T)| between zero and four (and more) times the modulus of 
(x)| as predicted by the chiral quark soliton model [26] yields the same results for {2Bqi^ussH^) within 1-cr uncer- 
tainty. This explains also why our results and those of Ref. [27] agree, as discussed in Sec. VII, in spite of the different 
models in particular for the d-quark transversity: hf{x) < in our approach [26] vs. hf{x) > in [27] from saturating 
the Soffer bound and choosing the positive sign. Sensitivity to the smaller (in the model [26] ) transversity antiquarks 
is far weaker. 

Notice that the Soffer bound severely constrains the region of r in Fig. 9. This means that the HERMES data 
[18, 19] leave hf{x) practically unconstrained. The Soffer bound provides basically all we know at present about the 
d-quark transversity, see Fig. 10. 



T.±(1/2X 



favoured 



unfavoured 



red , I 



-4 -2 2 4 r 

FIG. 9: (2_BGaussJi'i") as function of the param- 
eter r defined in Eq. (30). The Figure demon- 
strates that the SIDIS data on the proton Collins 
SSA [18, 19] are practically insensitive to hf{x). 
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The situation is different for the u-quark transversity. 
Varying ft," {x) within 30% of the chiral quark sohton model 
prediction (which is within the model accuracy) would re- 
sult in a 30% variation of both the favoured and unfavoured 
{2Bqi,ussH^) . Variations much larger than that would spoil 
the good agreement with the BELLE data in Sec. IV. Thus, 
the model prediction [26] for hf^x) seems to fit — within 
its theoretical uncertainty and within the uncertainties of 
our study — precisely in what is needed to explain simulta- 
neously the HERMES [18, 19] and BELLE [21] data with a 
compatible Collins fragmentation function. This gives rise 
to the picture for hi{x) shown in Fig. 10. 

To conclude the HERMES data [18, 19] on the proton 
Colhns SSA provide constraints for hi{x) but leave hf{x) 
practically unconstrained. This is not unexpected consider- 
ing the M-quark dominance in the proton. In order to learn 
more about hf{x) and eventually face the issue of a flavour 
separation data from different targets are required. The 
Collins SSA from a deuteron target is smaller than from a 
proton target and more difficult to measure as was shown 
by COMPASS [20]. It might be promising to explore a ^He 
target, see below Sec. X. 
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FIG. 10: xh\{x) vs. x in the kinematic range of the HER- 
MES experiment. The dashed/dotted lines are predictions 
for hi/hf from the chiral quark soliton model [26]. The 
shaded range for the ii-quark is bound from above by the 
Soffer bound and from below by the estimated theoretical 
uncertainty of the model result [26]. The shaded range for 
the d-quark is due to the Soffer bound. The Figure demon- 
strates roughly the emerging picture of the transversity dis- 
tribution. 



X. PREDICTIONS FOR ONGOING AND FUTURE EXPERIMENTS 



Further data on transverse target SSA in SIDIS is expected from COMPASS, HERMES and the JLab experiments. 
COMPASS plans to measure with a proton target, and HERMES with a deuteron target. Thus, the two experiments 
will supplement the results [17-20] in the respectively different kinematical regions. The best fit results to the data 
[17-20] in Figs. 3 and 7 show basically the predictions from our study for these experiments for charged pions. 

The Collins SSA for neutral pions is an interesting observable. Due to isospin symmetry this observable provides 
information which is already contained in the charged pions SSAs. However, being an independent experimental 
information it provides a valuable cross check for our understanding of the Collins function. If our present picture 
of the Collins function is correct, one may expect the 7r° Collins SSA to be rather small. Our prediction for the 7r° 
Collins effect on a proton target is shown in Fig. 11a. (On a deuteron target the effect is even smaller.) Preliminary 
HERMES data for the 7r° Collins SSA were shown in [17] and indicate an effect compatible with zero within large 
uncertainties. 

Rather interesting would be data on kaon Collins SSA. These would not only allow to test whether pions and kaons 
exhibit as Goldstone bosons of chiral symmetry breaking a similar Collins effect, see Eq. (28) and the subsequent 
discussion. The understanding of kaon Collins fragmentation is required to test the Schafer-Teryaev sum rule [73]. 

The Collins SSA for pions can also be measured in the CLAS and HALL-A experiments at JLab on different targets 
in a somehow different kinematics at ^ 2GeV^ and for 0.15 < a; < 0.5 which would provide further important 
information [75, 76]. For example, the HALL-A collaboration will measure with a transversely polarized "^He target 
which — after nuclear corrections — will provide data on the neutron Collins (and Sivers) effect. On the basis of the 
numbers in Eqs. (10, 11) we estimate for the pion Collins SSA on a "neutron target" the results shown in Figs, lib and 
11c as shaded regions. The estimates of the SSA could be even more optimistic because in the HALL-A experiment 
larger z = (0.4 — 0.6) will be probed than at HERMES, and the effect tends to increase at larger z, see Fig. 7. The 
displayed "error bars" in Fig. lib arc projections for 24 days of beam time from Ref. [76]. 

The TT^ Collins SSA on neutron is of similar magnitude as the tt^ Collins SSA on proton. However, the pion 
Collins SSAs on neutron are far more sensitive to hf{x). In order to illustrate this point we present in Figs, lib and 
11c (dotted lines) the result which follows from assuming that hf{x) has opposite sign to the chiral quark soliton 
model prediction [26] (in which hf{x) is negative). Clearly, on the basis of the HALL-A data it will be possible to 
discriminate the different scenarios, and to constrain hf{x) more strongly. HALL-A will also be able to measure the 
Sivers effect for kaons [76]. 

Before discussing the predictions for azimuthal asymmetries in e+e^ annihilation (via 7*) let us first establish how 
many different observables exist in the fragmentation into pions. There are free independent obscrvablcs, namely 

ylf(zi,Z2) - = ={^.^ + ^,^1 (31) 

A^{zi,Z2) = = ={^,^ + ^.^} (32) 



(33) 



11 



^si„»+*s)(x) on proton (a) AS![.<*+*s>(x) on neutron (b) A^![.<1'+*s)(x) on neutron (c) 




0.1 0.2 0.3 0.4 X 0.1 0.2 0.3 0.4 x 0.1 0.2 0.3 0.4 x 



FIG. 11: Ai"^'^+'^''^ vs. X as predicted on the basis of our study of the HERMES data in Sect. II for: (a) -k^ production from 
proton, which could be measured at HERMES [17] and CLAS [75]. (b) 7r~ production from a neutron (in practice '^He) target, 
as planned in the HALL-A experiment at JLab. The error bars estimated for the HALL-A experiment are taken from Ref. [76]. 
The dashed lines indicate tentatively the effects of a different scenario for hf{x), see text, (c) The same as (b) but for tv^ . 



Pj.(Zj, Zj) for 0.2 < z, < 0.3 (a) Pc'^l' ^2^ * ^1 * (b) ^c'^l' ^2) " * ^1 * ''c'^l' ^2* < < 1 (d) 




0.2 0.4 0.6 0.8 zj 0.2 0.4 0.6 0.8 Z2 0.2 0.4 0.6 0.8 Zj 0.2 0.4 0.6 0.8 



FIG. 12: The observable Pc{z\,Z2) as defined in Eqs. (35) and (36) for fixed 21-bins as function of Z2. The predictions are 
obtained on the basis of the fit result shown in Fig. 5, and could be tested in the BELLE experiment. 



The arrows indicate the opposite jets to which the pions tt^ belong, which carry the momentum fractions Zi of the 
fragmented quarks. (We discard the less convenient possibility of combining the events for charged and neutral pions.) 
These asymmetries are given by 

AY = l + cos(20i) /, {SiJffiJ^- + hH^fH^f + 2Hl^Hl^]/{H^^ ^ Dl) 
A[ = l + cos(20i) fe {^Hi^^H^f + ^H^fH^: + 2H:^HI^} / {H'^ ^ Dl} 

= 1 + cos(20i) /, {5(i?fr + H:f){Hi^: + n^f) + ah:^h:j/{h^ ^ on m 

where fg = sin^ 6'2/(l + cos^ 62) and = CGa_ussHi^^^^''"'{z) in our approach. 

The unlike- and like-sign asymmetries and were measured at BELLE (more precisely, the ratio A^ j A\ was 
measured) [21], sec Sect. III. The third asymmetry AS( can be accessed by observing a neutral and a charged pion 
as shown in the first line of Eq. (33). Actually, a linear combination of the observables in the first line of Eq. (33) 
contains the same information. Noteworthy, thanks to isospin symmetry the same information is contained in the 
asymmetry obtained by summing up all charged pion events — as indicated in the second line of Eq. (33). Thus, this 
alternative way allows to access the information content of possibly more easily — considering the more difficult 
7r° detection. 

To get rid of hard gluon and detector effects, sec Sect. III. one may consider the double ratio 



Ac 



l + cos(20o) ^'c(^i,^2) , (35) 
where Pc{zi, Z2) is given to leading order in the analyzing power by 

Pc{zi,Z2) = fe 



(36) 



in the notation introduced in the context of Eq. (34). 

The predictions for Pc(-zi, ^2) made on the basis of our study of the BELLE data [21] in Sect. Ill are shown in Fig. 12. 
We observe a qualitatively similar picture to that in Fig. 6 for the observable Pi{zi, Z2) defined in Eqs. (16, 17). Thus, 
Pc{zi,Z2) could be seen in the BELLE experiment in a similarly clear way as Pi(zi,Z2)- We stress that Pc{zi,Z2) 
contains independent experimental information which would provide valuable constraints allowing to better pin down 
the Collins function and to discriminate between H^^'" and iJj*-""^ using less constrained Ansatze than e.g. (18). 
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XI. CONCLUSIONS 

We presented a study of data on azimuthal asymmetries due to the Collins effect in SIDIS and in e+e^ annihilation. 
Our investigation bears considerable theoretical uncertainties due to the neglect of soft factors, Sudakov effects, scale 
dependence and choice of models for transverse parton momenta and the transversity distribution. We tried to take 
into account these uncertainties and to minimize them wherever possible. For example in order to confront results from 
experiments performed at different scales we preferred to compare appropriately defined ratios of spin-dependent to 
spin-independent quantities which might be expected to be less scale dependent than the respective absolute numbers. 

We observed that — within the uncertainties of our study — the SIDIS data from HERMES [17-19] and COMPASS 
[20] on the Collins SSA from proton and deuteron targets are in agreement with each other and with BELLE data 
[21] on the azimuthal asymmetry in e'^e^ annihilation, and seem compatible with the indications from DELPHI [16] 
for an azimuthal asymmetry in e~'"e~ annihilation at the Z'^-peak. 

The emerging picture for the so far unknown functions is as follows. 

• The Collins fragmentation functions show the behaviour H^^^^^^ oc zDi(z) at small z as found here and in [27]. 
Within our model for transverse momenta this relation remains valid also at large z, but not in the approach 
of [27] . The favoured and unfavoured Collins fragmentation functions are of comparable magnitude but have 
opposite sign. SIDIS data require the favoured Collins fragmentation to be positive in the convention [77]. 

• The u-quark transversity distribution is positive and close (within 30%) to saturating the Soffer bound. In 
contrast, hf{x) and antiquark transversities are entirely unconstrained by present data. At the present stage 
the chiral quark soliton model [26] provides — within its theoretical uncertainty — a useful estimate for hi(x). 

Future data from SIDIS and e'^e~ colliders will help to refine and improve this first picture. Useful for that will be 
the confrontation of these data with the estimates for observables in SIDIS and e+e^ annihilation we presented which 
can be measured in running and planned experiments at HERMES, COMPASS, CLASS, HALL-A and BELLE. 

Particularly helpful would be data on double transverse spin asymmetries in Drell-Yan Icpton pair production from 
RHIC [78] and the planed PAX experiment at GSI [79] which could provide direct information on hf{x). In particular, 
the PAX experiment — making benefit from a polarized antiproton beam — could provide valuable constraints on the 
M-quark transversity distribution [80]. Only the combined analysis of SIDIS, e~^e~ and DrcU-Yan data will provide 
the full picture of the Collins fragmentation and transversity distribution function. 
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